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ABSTRACT 

We investigate the physical properties, geometry and dynamics of the massive cluster 
merger MACS JO 140.0-0555 (z=0.451) using X-ray and optical diagnostics. Featuring two 
galaxy overdensities separated by about 250 kpc in projection on the sky, and a single peak in 
the X-ray surface brightness distribution located between them, MACS JO 140.0-0555 shows 
the tell-tale X-ray/optical morphology of a binary, post-collision merger. Our spectral analysis 
of the X-ray emission, as measured by our Chandra ACIS-I observation of the system, finds 
the intra-cluster medium to be close to isothermal (~8.5 keV) with no clear signs of cool 
cores or shock fronts. Spectroscopic follow-up of galaxies in the field of MACS J0140.0- 
0555 yields a velocity dispersion of 875^q km s _1 (n z = 66) and no significant evidence of 
bimodality or substructure along the line of sight. In addition, the difference in radial velocity 
between the brightest cluster galaxies of the two sub-clusters of 144 ± 25 km s _1 is small 
compared to typical collision velocities of several 1000 km s _1 . A strongly lensed background 
galaxy at z=0.873 (which features variable X-ray emission from an active nucleus) provides 
the main constraint on the mass distribution of the system. We measure M(<75kpc) = (5.6± 
0.5) x 10 13 M for the north-western cluster component and a much less certain estimate of 
(1.5 — 3) x 10 13 M for the south-eastern subcluster. These values are in good agreement with 
our X-ray mass estimates which yield a total mass of MACS J0140.0-0555 of M(<r 5 oo) ~ 
(6.8-9.1) x 10 14 M . 

Although all optical and X-ray properties of MACS JO 140. 0-05 5 5 are consistent with a 
well advanced head-on merger proceeding along an axis close to the plane of the sky, the 
degeneracy between Hubble flow and peculiar velocity prevents us from obtaining a quantita- 
tive constraint on the inclination angle of the merger axis. The lack of pronounced substruc- 
ture in the cluster gas distribution and the proximity of the two optical cluster cores both in 
projection and in radial velocity suggests that the merger is observed well after the primary 
collision and, possibly, after turnaround. A weak-lensing analysis of the mass distribution in 
MACS JO 140.0-0555 has the potential of yielding constraints on the self-interaction cross sec- 
tion of dark matter similar to those obtained for MACS J 0025.4—1222, although the smaller 
projected separation of the two cluster components makes the measurement more challenging. 

Key words: Galaxy clusters: individual: MACS J0140.0-0555 - 



1 INTRODUCTION 

As the largest gravitationally bound systems in the universe, clus- 
ters of galaxies harbour copious amounts of both luminous and dark 
matter. Located at the densest regions of the Cosmic Web ( Bo nd et| 



|* Some of the data presented herein were obtained at the W.M. Keck Ob-| 
servatory, which is operated as a scientific partnership among the California 
Institute of Technology, the University of California and the National Aero- 
nautics and Space Administration. The Observatory was made possible by 
the generous financial support of the W.M. Keck Foundation. 



|al.|1996| ), clusters grow through semi-continuous accretion of mat- 
ter along large-scale filaments as well as through discrete mergers 
with nearby clusters. The dynamics of both of these processes are 
well explored by numerical simulations (e.g.,| Roettiger et al. 1996 
Colberg et al.|1999| [2005] |Ritchie & Tliomas|2002| |CeccarellTet| 
al.|2011| >. Detailed studies of cluster mergers thus provide a unique 
opportunity for studies of the physics of mass assembly at the nodes 
of the Cosmic Web, as well as of the dynamical behaviour of dark 
matter. 

A crucial and obvious prerequisite for such studies is the iden- 
tification of massive cluster mergers, either active or very recent. A 
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shortlist of extreme, but geometrically simple, linear mergers was 
recently compiled by Mann & Ebeling (2012 ) by examining the X- 
ray and optical morphology of a statistical sample of very X-ray 
luminous clusters at z > 0.15. By virtue of being binary, head- 
on mergers, these systems should be particularly well suited for 
quantitative interpretation of the immediate observables. The latter 
include a simple, yet powerful diagnostic that can be used to in- 
vestigate the nature and properties of dark matter, namely the spa- 
tial separation of collisional and non-collisional cluster components 
during cluster mergers. If dark matter is indeed collisionless (e.g., 
Gne din & Ostr iker 2001, Natarajan et al. 2002), cluster mergers 
should result in a clear segregation between the collisionless dark 
matter and the viscous intra-cluster medium (ICM) which is subject 
to shock heating and compression. 

Extreme cases of a segregation of this kind were detected 
in 1E0657-56 0=0.296; the Bullet Cluster; |Clowe et al.l|2004l 
12006] ), and MACS J0025 .4-1222 (z=0.586; [Bradac et al.|2008) . 
Although alternative explanations are still being explored (e.g., 
modified Newtonian dynamics and tensor-vector-scalar gravity; 
|Angus et al.|2007||FeTx et al.|20Q8) , the observed offsets between 
the centroid of the X-ray emission and the distribution of the to- 
tal gravitational mass as reconstructed from weak-lensing observa- 
tions provides strong and direct empirical evidence of the existence 
of dark matter, and allows upper limits to be placed on the dark 



matter self-interaction cross-section, a/m ([Markevitch et al. 2004 



|Randall et al."l2008) . Of the various methods available to place ob 
servational limits on a/m from cluster mergers, direct measure- 
ments of the mass surface density distribution in a post-collision, 
head-on, equal-mass merger are perhaps the most powerful as they 
can be obtained straightforwardly for any massive merger, as long 
as the merging components are well separated on the sky (Marke- 
|vitch et al.|2004||Bradac et al.|2008) . 

We here use X-ray and optical observations of 
MACSJ0140.0-0555, one system from the list of |Mann &| 
Ebeling (2012), in order to unravel the dynamics and geometry 
of the merger, and to establish its usefulness for measurements of 
a/m. 

1.1 MACSJ0140.0-0555 

MACSJ0140.0-055 5 Q=0.451) is one of several c lusters from 
the MACS sample ( [Ebeling et al"][200T1 [2007] [20T0] > with an X- 
ray/optical morphology similar to that of MACS J0025. 4-1 222. 
In the optical, MACS J0140.0-0555 clearly comprises two sub- 
clusters separated by about 250 kpc in projection on the sky. At 
X-ray wavelengths, a short (10 ks) Chandra observation revealed 
extended emission from the gaseous ICM elongated along the axis 
connecting the brightest cluster galaxies (BCGs) of the two sub- 
clusters and featuring a single peak between the BCGs (Fig. [TJ. 
Consequently, the system was selected by Mann & Ebeling (2012 ) 
as one of eleven likely post-collision, binary, head-on mergers. A 
snapshot image of the system (see inset in Fig. [TJ, taken with the 
Advanced Camera for Surveys (ACS) aboard the Hubble Space 
Telescope (HST) shows a putative gravitational arc close to the peak 
of the extended ICM emission. If confirmed, this arc would yield 
immediate constraints on the mass of one, and possibly of both of 
the two apparent subclusters. 

In order to confirm this tentative picture and to obtain a quan- 
titative description of the merger history and physical properties of 
MACS J0140.0-0555, we conducted additional observations with 
Chandra and HST, as well as groundbased galaxy spectroscopy 
with Keck telescopes. In this paper, we present the results of our 
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Figure 1. Iso-intensity contours of the adaptively smoothed (Ebeli ng et al.| 
2006 ) X-ray emission from MACS J0140.0-0555 as observed with Chan- 
dra (ACIS-I; 10 ks exposure) overlaid on an optical image (V, R, I) obtained 
with the UH2.2m telescope. Contour levels are spaced logarithmically. The 
inset shows a putative gravitational arc near the peak of the X-ray emis- 
sion, detected by an HST snapshot image taken for program GO- 10491 (PI 
Ebeling). 



comprehensive analysis of all data collected. Specifically, we per- 
form spatial and spectral analyses of the Chandra X-ray data to 
quantify the distribution, mass, and temperature of the ICM. In 
addition, we assess the distribution and dynamics of the cluster 
galaxy population based on groundbased spectroscopy. The latter is 
also used to obtain redshifts of putative gravitational arcs that may 
yield strong-lensing constraints on the dark matter distribution in 
MACS J0140.0-0555. The paper is structured as follows. Section 
2 describes the observations and data reduction. Section 3 presents 
the results from of our X-ray and optical analyses. In Section 4, 
we discuss these results and their implications for the nature of the 
merger. A summary is given in Section 5. 

Throughout this paper, we assume the concordance ACDM 
cosmology with H = 70 km s _1 Mpc -1 , Qm = 0.3, and 

= 0.7, for which \" corresponds to 5.77 kpc at the cluster 
redshift of z=0.451. In addition, abundances are measured using 
the abundance table supplied by A nders, E. & Grevesse, N.| fl989 ). 
All errors are la unless stated otherwise. 



2 OBSERVATION AND DATA REDUCTION 
2.1 X-ray data 

MACS J0140.0-0555 was observed with the A dvanced CCD Imag- 
ing Spectrometer (ACIS,[Garmire et al. 2003 ) aboard the Chandra 
X-ray observatory in 2004 June (ObsID 5013) and 2010 Septem- 
ber (ObsID 12243) for a total exposure time of 10.3 ks and 19.6 ks 
respectively in VFaint mode. The data were reduced following 
standard Chandra data reduction procedure^ using CIAO 4.3 and 
CALDB 4.4.2 ( Frusc iorie et al.|2006| ). No filtering was necessary as 
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no obvious flares occurred during either observation. For the spatial 
analysis described in detail below, we merged the two event files af- 
ter proper alignment. The spectral fitting, however, was performed 
separately (but simultaneously) on the two event files. 

In order to avoid instrumental artifacts caused by position- 
dependent charge transfer inefficiency, we did not measure the 
background from an off-source region of our actual datasets, but 
instead used "blank-sky" backgrounds drawn from the calibration 
database. We ensured that the backgrounds from the "blank-sky" 
are at the same levels as the observations by comparing photon 
count rates at high energy (10.1-1 1 keV) where emission predom- 
inately comes from background. The background datasets for the 
two observations were treated as the science datasets, i.e. they were 
merged for the spatial analysis, but dealt with separately during the 
spectral fitting. 

In all fits to the data, spatial or spectral, point sources were 
excluded. 



2.2 Optical data 

2.2.7 Hubble Space Telescope imaging 

MACSJ0140.0-0555 was observed in the F606W filter with the 
Advanced Camera for Surveys (ACS, |Giavalisco et al.| [2004) 
aboard the Hubble Space Telescope on October 3, 2005 for 1,200 
seconds as part of Snapshot programme GO-10491 (PI: Ebeling). 
The data were corrected for CTI effects using the pixel-based al- 
gorithm developed by Andreson & Bedlin (2010) for this pur- 
pose, then multi-drizzled (Koekemoer et al. 2002) and resampled 
( [Rhodes et al.| |2007) into 0.03"pixels. A second observation was 
performed on September 12, 2011 with the Wide Field Camera 3 
(WFC3) in the F110W and F140W filters both for 706 seconds as 
part of Snapshot programme GO- 12 166 (PI Ebeling). 



2.2.2 Groundbased spectroscopy 

Follow-up spectroscopy of likely cluster galaxies and potential 
strong-lensing features in MACS J0140.0-0555 were performed 
with the Low Resolution Imaging Spectrograph (LRIS, Oke et al.] 
[T995l|McCarthy et al.|1998||Rockosi et al.|2010) on the 10m Keck- 
I telescope on Mauna Kea on December 10, 2010. Using the 600 
line/mm grating blazed at 7500A on the red arm of the spectro- 
graph, the 6800A dichroic as a beam splitter, and the 300/5000 
grism on the blue arm, a multi-object spectroscopy mask with 
38 science slits was observed for 3 x 1900 seconds. Standard 
data reduction procedures (bias subtraction, flat fielding, sky sub- 
traction, wavelength calibration) were applied before the extrac- 
tion of one-dimensional spectra for each target object. Twenty 
eight objects yielded useful spectra for further analyses. Additional 
spectroscopic follow-up observations were performed on Septem- 
ber 11, 2011 with the DEep Imaging Multi- Object Spectrograph 
(DEIMOS; Faber et al. 2003) on Keck-II telescope, where we used 
the GG455 order-blocking filter and the Zerodur 600 line/mm grat- 
ing centred at 6300A. Four slit masks containing a total of 60 
science slits were designed for this observation and observed for 
3 x 600s; 44 of them yielded useful spectra after standard data 
reduction procedures were applied. Both data reductions are per- 
formed using our pipeline adapted from the Deep Extragalactic 
Evolutionary Probe 2 Survey data reduction pipeline ( Davis ^TaT] 
|2003] ). 




0.1 1.0 

Radius (arcmin) 



Figure 2. Photon flux as a function of radius measured from the peak of the 
diffuse X-ray emission. The red line indicates the best- fit 1-D /3-model (see 
Section 3.1) with an additional background term. The dashed line marks the 
best-fit background level. 



Table 1. Best fit parameters of 2-D (3 -model 



r e 1 (5 

2033 i 14.0 k P c 0.24 ±0.02 27.2° ± 2.4 0.63 ± 0.03 



e = 1 — b/a; a and b are the major and minor axis, respectively. 

3 DATA ANALYSIS AND RESULTS 
3.1 X-ray surface brightness 

We quantify the X-ray surface brightness distribution (see Fig. [TJ 
by fitting the data with the "isothermal King" model (also known as 
/3-model; Cavaliere & Fusco-Femi ano| 1 976} . In the simple spher- 
ical case, the gas density profile p(r) and the projected surface 
brightness profile S(r) can then be described by 

p (r)= P0 [l + (I-f ] -¥ 

and 

5(f) = 5 [l + (-) 2 ]- 3/3+0 - 5 5 
ro 

where r is radius r in projection on the sky, p is the central gas 
density, So is the central surface brightness, f3 is the power index, 
and ro is the characteristic core radius. 

The observed one-dimensional surface brightness profile 
(Fig. [2} is computed as a function of radius by measuring the sur- 
face brightness in circular annuli centred on the peak of the X-ray 
surface brightness. Fitting a 1-D /3-model yields f3 — 0.761q'°q 
and ro = 209^24 kpc. The goodness of fit as given by the reduced 
X 2 statistic is 1.06. 

Since the obvious ellipticity of the X-ray emission can not 
be well described by the circular 1-D profile above, we use, for 
our two-dimensional model, a more general, elliptical /3-model 
(beta2d in Sherpa), which adds two more fit parameters, namely 
the ellipticity e and the position angle 0. The resulting best-fit pa- 
rameters are given in Table [T] and the best-fit model and residu- 
als are shown in Fig. [3] Because of the poor photon statistics at 
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Figure 3. X-ray surface brightness distribution (left), 2D /5-model of the distribution (middle ), and residuals (r ight). A constant background is included in the 
fit. The data and residuals are binned into 16" pixel. Contours show the adaptively smoothed ( Ebeling et al. 2006 > X-ray emission. Note the obvious excess in 
the residuals located south-east to the peak of the X-ray surface brightness. 



large cluster-centric radii, the Cash statistic ( Cash 1979) is adopted; 
therefore a formal goodness of fit is not available. However, a sim- 
ple visual inspection of Fig. [3] proves instructive. While the over- 
all elongation of the emission is well captured by our elliptical f3- 
model, the deviation from concentricity in the cluster core leads to 
pronounced residuals when the model is subtracted from the data 
(right panel of Fig. [3J. Comparing the residuals to the raw data (left 
panel of Fig. [3} and to the X-ray /optical overlay shown in Fig. [T] 
we find these residuals to be caused by the ridge-like structure in 
the X-ray surface brightness that spatially coincides with the core 
of the south-eastern (SE) sub-cluster. 



3.2 Global X-ray properties 

Global X-ray properties are measured within risoo (0.68 Mpc or 
118"f] which is the radius at which the signal-to-noise ratio (S/N) 
of the global cluster emission approximately reaches its maximum, 
such that extending the measurement to yet larger radii would pre- 
dominantly add background (Fig. [2]). Inside risoo, we measure a 
S/N of 69 for a total of 5445 net X-ray photons in the 0.3-6 keV 
band. For the spectral fit, we adopt the plasma model of |Mewe| 
et al.|(|1985| combined w ith the photoelectric absorption model of 
Morrison & M cCammon| { |1983| ) (XSMEKAL and XSWABS, respec- 
tively, in Sherpa); all fits are performed simultaneously on the two 
observational datasets available for MACS J0140.0-0555. We find 



n H = (4.5 ± 2.0) x 1 20 cm 2 , consistent with the Galactic value 



of 2.85 x 10^ u cm" 2 ( [Dickey & Lockman|l990) ; a global temper- 
ature of 7.7^0 keV; and a metal abundance of Z = 0.35 ± 0.12, 
consistent with results obtained for other clusters at similar redshift 
(Arn aud et al.|1992) . The goodness of fit of a single-temperature 
model characterized by a reduced \ 2 of 104 confirms the reason- 
ableness of the isothermal assumption. Adding a second tempera- 
ture component to the model does not yield a better fit. In order 
to obtain better constraints on the ICM temperature, we freeze the 
absorption at the Galactic value and also set the poorly constrained 



2 ?"1500 is the radius at which the mass overdensity of the cluster exceeds 
1500 times the critical density. Like the similarly defined radii r2oo (1.9 
Mpc, 324") and rsoo (1-2 Mpc, 205"), risoo is computed as a function of 
kT and z ( |Ma et al.|2008) . 



global abundance to Z — 0.3. This yields a global gas temperature 
of8.5t^keV. 

For this final set of model parameters, the total rest-frame X- 
ray luminosity of the cluster within risoo is measured to be 1.0 x 
10 45 erg s"\ 4.9 x 10 44 erg s~\ 9.6 x 10 44 erg s~\ and 1.6 x 
10 45 erg s" 1 (0.3-6 keV, 0. 1-2.4 keV, 2-10keV and bolometric, 
respectively.). Extrapolating these figures to rsoo using the 2D /3- 
model, we find values for Lx(< rsoo) of 1.3 x 10 45 erg s _1 , 
6.5 x 10 44 erg s"\ 1.3 x 10 45 erg s~\ and 2.1 x 10 45 erg s" 1 
in the same four passbands. 

The normalization factor of the global spectral fit allows us to 
derive an estimate of the total X-ray gas mass within the region, 
in as much as it yields the central electron density that permits the 
conversion of the projected X-ray surface brightness into a phys- 
ical, three-dimensional gas density. The spectral normalization of 
the XSMEKAL model is given by 

10 " 14 f\ Ue n Uh \j\r 
norm = p. / \dV. 

^[%t\ 2 (l + z) 2 J cm3 cm3 

We calculate the right-hand side of this expression by numeri- 
cally integrating the hydrogen number-density profile. The gas is 
assumed to be fully ionized, yielding the electron number den- 
sity, n e , of 1.2 times the hydrogen number density. The hydrogen 
and electron density profiles are de-projected from the X-ray sur- 
face brightness fit, i.e. 2D elliptical /3-mode|^]as described in Sec- 
tion |TT] Taking into account the known steepening of the X-ray 
surface brightness profile at large radii where our data do not con- 
strain the mass density model, we follow Vikhli nm et al.| (20Q6) 
and adopt a steeper slope of 0.60 at r > 0.3 x r2oo and a den- 
sity cutoff at r > rsoo- Without these modifications, the integral 
over the best-fit /3 model does not converge. Using this modified 
model, we find values for the central electron number density, n e , 
of 8.9 x 10 -3 cm -3 and for the total X-ray gas mass within risoo of 
(6.55 ±0.01) x 10 13 M©. The quoted uncertainties are the nominal 
la errors propagated from the normalization in the spectral fitting. 
Note, however, that the real uncertainties are likely to be dominated 
by systematic errors, since we assume that (1) the gas is isothermal, 
(2) the gas body is a prolate spheroid with its major axis lying in 



in doing so, we assume a prolate geometry 
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spectral fits are performed in each region as described in Section 
3.2| Of the model parameters, only the gas temperature and the 
normalization are fit, while the redshift is fixed at the cluster red- 
shift of 0.451, the absorption term is set to t he Galactic value of 
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Figure 4. Profiles of the cylindrically enclosed (2D) and spherically en- 
closed (3D) gas mass derived for different geometries and X-ray surface- 
brightness profile slopes. For each set of curves, the upper two assume a 
slope of 0.6 for the X-ray surface brightness profile at large radii, and the 
lower two assume 0.92; the uncertainties caused by the adopted geome- 
try are negligible in comparison as an oblate geometry results in an only 
slightly higher mass than a prolate one. Details are provided in Section[T2] 



the plane of the sky, and (3) the density profile is described by a 
modified /3-model. While the assumption of isothermality on large 
scales is supported by our analysis of the ICM temperature distri- 
bution (Section [33] ), the 3-D geometry of the X-ray gas is poorly 
constrained, and the /3-model provides a poor fit to the central re- 
gion of the X-ray surface brightness distribution (Figs. [2] [3]). 

To understand the systematic errors caused by the unknown 
geometry and uncertain profile at large radii, we repeat the calcu- 
lation for an oblate spheroid (with a major axis that lies again in 
the plane of the sky) and also vary the slope of the X-ray surface 
brightness at r > 0.3 x r2oo from 0.60 to 0.92. The results are 
shown in Fig|4] Two different types of gas mass as a function of r 
are shown; the gas mass enclosed by an infinitely long cylinder of 
radius r pointed toward the observer and those enclosed by a sphere 
of radius r (the 2-dimensional and 3-dimensional gas mass, respec- 
tively). The 2D and 3D gas mass profiles illustrate that systematic 
effects indeed dominate at radii greater than risoo, especially for 
the 2D gas mass. We find M(< n 5 oo) ~ (5.7 - 6.9) x 10 13 M 
(2D) and M(< n 5 oo) ~ (4.5 - 4.7) x 1O 13 M (3D). At r 50 o 
systematic errors increase dramatically, leading to a more uncer- 
tain measurement of M gas (< ^soo) = (7.8 - 10.5) x 10 13 M 
(2D and 3D are alike due to density cutoff at r > rsoo). The mass 
ranges quoted here are the maximum and minimum enclosed gas 
mass of the four different types of geometries. 

The total gas mass within rsoo derived above implies a total 
gravitational mass for MACS J0140.0-0555 of M to t(< rsoo' 

9.1) xlO 14 M , assuming a baryon fraction of 0. 1 15 (Mantz 



(6 

|et al.|2010|| Allen et' al. 2008 ) and neglecting the mass in stars 



3.3 X-ray temperature map 

In recognition of the limited photon statistics, we adopt the two- 
level binning method described in |Ma et al. (2009 ) for our analysis 
of ICM temperature variations. First, sub-regions are defined in the 
count image (0.3-6 keV) using the "cumulative binning" algorithm 
(contbin; Sander 2006) for a given S/N threshold. Subsequently, 



2.85 x 10 20 cm" 2 ( [Dickey & Lockman 



1990 ), and the metal abun- 



dance is frozen at a value of 0.3, typical of clusters ( [Arnaud et al.| 
[T992| [Balestra et al.|2007] ). A S/N criterion of 15, corresponding 
to typically ~ 270 net photons in each region, yields a first high- 
resolution, low -fidelity temperature map. We then manually com- 
bine adjacent regions of similar temperature, thereby creating areas 
comprising up to 1,300 net photons, and perform the spectral fits 
again to arrive at a low-resolution, high-fidelity temperature map. 

The top panels of Fig.|5]show the low-fidelity temperature map 
and the temperatures measured in each individual region. Within 
the large measurement errors, the temperatures of all regions are 
consistent. However, careful inspection of the map suggests that 
the regions to the west and north-west (labeled 10, 7, 6) have con- 
sistently higher temperatures than the regions to the east and south- 
east (i.e., 0, 1,2, 3). We therefore rebin these regions to increase 
the photon statistics. Examining the rebinned map, we find the dif- 
ference in the ICM temperatures of the merged regions to be more 
pronounced but still statistically insignificant (i.e. regions labeled 
A, B). 

As an alternative probe of the isothermality of the ICM we 
also measure the gas temperature in concentric annuli, i.e., as a 
function of radius. Requiring approximately 1100 net photons per 
annulus, we arrive at the temperature map and profile shown in 
the bottom panels of Fig. [5] Again we find no significant devia- 
tions from the average gas temperature of approximately 8.5 keV, 
but note that variations at the level suggested by the results shown 
in the middle panel of Fig. [5] may still be present but undetectable 
with the current photon statistics of our data. 

3.4 Redshifts 

The reduced 1-D spectra of all galaxies and gravitational arc can- 
didates targeted in our spectroscopical follow-up observations are 
analyzed with the IDL routine SpecPro developed by Masters 
& Capak (2011). Spectroscopic redshifts are determined by cross- 
correlation against spectral templates. Given that most of the galax- 
ies are red ellipticals, we use the VVDS E0 template, except for 
five objects with strong narrow or broad emission lines for which 
the VVDS Starburst and SDSS Quasar templates, respectively, are 
adopted. The results are tabulated in Table [2] 

Our redshift measurements allow us to identify 66 cluster 
members with redshifts between 0.423 and 0.461, two foreground 
objects, and four background objects, one of which is the gravita- 
tional arc (Al) discussed in more detail in Section [33] Using the 
biweight and bootstrapping estimator from the ROSTAT statistics 
package of Bee rset al.| ([T990 ) and converting redshifts to heliocen- 
tric radial velocities (including relativistic corrections), we calcu- 
late a mean cluster redshift of 0.4465±£JJj£g and a radial velocity 
dispersion of 875^q km s _1 . A histogram of the redshift dis- 
tribution is shown in Fig. [6] Within the limits set by the number 
statistics, the distribution is well described by a Gaussian: a two- 
sided Kolmogorov-Smirnov test yields a probability of 88% for the 
observed distribution to be drawn from a normally distributed par- 
ent distribution. Visual inspection of the distribution of redshifts on 
the sky (Fig. [7} shows no obvious correlation between the positions 
and the redshifts of the cluster members. 

We measure the redshifts of the two BCGs as 0.4502±0.0001 
for the north-western BCG (BCG1) and 0.4510 ± 0.0001 for the 



© 0000 RAS, MNRAS 000, 000-000 



6 Ho, Ebeling & Richard 



-05° 57 



53 




25 



54- 



55 



56 



-05° 57'- 



-05° 57 



05 s 40 m 00 s 01 n 39 m 55 s 

Right Ascension (2000) 



6.0 13.0 
Temperature (keV) 




05 s 40 m 00 s 01 h 39 m 55 s 

Right Ascension (2000) 




> 

CD 



20 



15 



10 



5 
1.2 
1.0 

CM 

* 0.8 
0.6 



25 



_ 20 
> 

CD 
CD 

=5 15 

CO 
CD 



10 



5 
1.2 
1.0 

CM 

* 0.8 
0.6 



25 



_ 20 
> 

CD 



5 15 



10 



5 
1.4 
1.2 
w x 1.0 
0.8 
0.6 



4 6 8 
Region number 



10 



12 



AB23456789 10 11 12 
Region number 



05 s 40 m 00 s 01 n 39 m 55 s 

Right Ascension (2000) 



-1 



1 

Region number 



Figure 5. Left panels: temperature of the gaseous ICM estimated by fitting plasma and photoelectric absorption models to the X-ray emission in different 
regions (details are provided in Section [33) . The temperature is color coded as indicated by the color bars. Blue contours show the adaptively smoothed X-ray 
emission. Right panels: the best-fit temperature and reduced x 2 value for each region labeled in the left panels. 
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Figure 6. Distribution of spectroscopic redshifts of galaxies in 
MACS J0140.0-0555 as determined from our LRIS and DEIMOS obser- 
vations. The errors and upper limits shown assume the number of sources 
in each bin follows poisson statistics. The dashed line describes a Gaus- 
sian distribution with a dispersion of 875 km s — 1 centred on the systemic 
cluster redshift of z = 0.4465. 



south-eastern BCG (BCG2). This corresponds to a difference in 
radial velocity of 144 ± 25 km s _1 . 



parametric models created with the Lenstool software ( Julio et aT] 

|2007l >. 

In our model of the mass distribution, we include both smooth 
cluster- scale components as well as galaxy- scale components. For 
the latter, each cluster galaxy is assigned a small-scale isothermal 



elliptical halo scaled according to its luminosity (see Richard et 
|al.|2010| for more details on the modeling process). For the cluster 
scale mass haloes, we assume either one pseudo-isothermal compo- 
nent centered on the north-western BCG, or two components cen- 
tered on either BCG. Although we cannot precisely constrain the 
mass ratio between the main and the secondary cluster-scale com- 
ponents, we find that the large radius of curvature of Al is better 
reproduced by models with two haloes, with the second, south- 
eastern halo containing up to 50% of the total mass of the main 
halo. The geometry of the critical lines (Fig. [TJ is also compatible 
with other strong-lensing features in MACS JO 140.0-0555, found 
to be single images of galaxies with spectroscopic redshifts of z~l. 

While these simple models provide no tight constraint on the 
mass of the south-eastern cluster component, they yield a robust 
measurement of the mass of the north-western main cluster. We find 
M(<75kpc) = (5.6 ±0.5) x 10 13 M for the projected mass out 
to the radial distance marked by the location of Al. This value is 
free of systematic uncertainties as long as the mass distribution of 
the north-western cluster is centred on the location of BCG1 within 
1", an assumption that is supported by strong-lensing analyses of 
many other clusters (Ric hard et al.|2010"| . 



3.5 Central gravitational arc 

The presence of a potential strong-lensing feature between the two 
BCGs, labeled Al in Table[2] raises the possibility of obtaining di- 
rect constraints on the mass distribution of MACS J0140.0-0555. 
In the optical (Fig. [7}, Al appears to be a nearly straight gravita- 
tional arc at an approximately tangential orientation with respect 
to the two BCGs. LRIS spectra of Al are shown in Fig. [8] with 
identified lines labeled. Strong emission lines in both the blue and 
red side of the spectrum allow an unambiguous redshift measure- 
ment of z=0.8729, placing this object well in the background of 
MACS JO 140. 0-05 5 5 and confirming its nature as a gravitational 
arc. The broad line profiles and the unresolved, variable X-ray 
emission suggest that Al is a Seyfert I active galactic nucleus. 

In the X-ray regime, Al coincides with a highly time- variable 
point source. We measure 0.3-6 keV fluxes of (4.9l4 3 9 5 ) x 
10 -16 erg s _1 cm -2 from the first Chandra observation in 2004 
(too little for the source to be visible in Fig.[TJ, and (1.1 ± 0.3) x 
10 -14 erg s _1 cm -2 from the second one in 2010. However, the 
implied variability of at least a factor of 4 (but conceivably much 
more) over a time of several years still needs to be adjusted by the 
gravitational amplification by the cluster at the location of the arc. 
Our mass model (see below) estimates this magnification factor to 
be about 10 =b 2 which reduces the variability to a modest level 
com monly observed in AGN (e.g.,[M arkowitz, Edelson & Vaughan 
|2003l|Markowitz & Edelson|2004[|Papadakis et al.|2008| >. 

3.6 Mass model from strong gravitational lensing 

Strong-lensing constraints are provided by the location, geometry, 
and orientation of Al, all of which identify it as a highly magni- 
fied tangential arc, either very close to or crossing the critical line. 
We attempt to reproduce the observed constraints using a series of 



4 DISCUSSION 

In the following we attempt to address three questions: a) what is 
the inclination of the merger axis relative to the plane of the sky, 
and can we constrain the impact parameter of the collision, b) what 
is the likely total mass of MACS J0140.0-0555 and where does the 
system fall on the scaling relations regarding mass, X-ray tempera- 
ture, and X-ray luminosity, and c) can this cluster be used to tighten 
the existing constraints on the self-interaction cross section of dark 
matter? 

4.1 Merger dynamics and geometry 

Our analysis of the imaging and spectroscopic evidence at both op- 
tical and X-ray wavelengths confirms our original hypothesis that 
MACS J0140.0-0555 is an active merger of two sub-clusters. A sin- 
gle emission peak is observed in the X-ray regime, located between 
the two BCGs, and no evidence of a surviving cool core is appar- 
ent in either the X-ray surface brightness or the ICM temperature 
maps. The presence of a gravitational arc with a spectroscopic red- 
shift of z = 0.8729 yields strong-lensing constraints on the mass 
distribution that favour a two-component mass model. The relative 
masses of the two subclusters are, however, poorly constrained but 
likely to fall in the range of 1:4 to 1:2. Spectroscopic redshifts of 
66 cluster members show no evidence of a bimodal distribution in- 
dicative of a line-of- sight merger and also no obvious correlation 
between radial velocity and position on the sky. Our tentative con- 
clusion that the merger axis lies close to the plane of the sky will 
be examined more closely in the remainder of this section. 

We first note that the relative radial velocity of the two BCGs 
of 144 zb 25 km s _1 is very small compared to the terminal ve- 
locity of BCGs in massive clusters undergoing a head-on collision. 
Markevitch et al. (2004) estimate the merger velocity for the Bul- 
let Cluster to be 4500 km s _1 by calculating the free-fall velocity 
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Table 2. Spectroscopic galaxy redshifts in the field of MACS J0140.0-0555. 



R.A. Dec. z R.A. Dec. z 



LRIS 



01 


39 


49.17 


-05 


52 


27.78 


0.4508 


01 


39 


57.34 


-05 


56 


17.42 


0.4383 


01 


39 


49.26 


-05 


52 


4.20 


0.4508 


01 


39 


57.35 


-05 


54 


34.73 


0.4368 


01 


39 


49.88 


-05 


52 


50.27 


0.4436 


01 


39 


57.44 


-05 


55 


1.12 


0.4484 


01 


39 


50.40 


-05 


52 


20.33 


0.3705 


01 


39 


58.03 


-05 


54 


39.57 


0.4494 


01 


39 


50.42 


-05 


55 


33.95 


0.4430 


01 


39 


58.21 


-05 


54 


57.98 


0.4610 


01 


39 


50.70 


-05 


54 


31.47 


0.4563 


01 


39 


58.38 


-05 


55 


8.30 


0.4577 


01 


39 


52.62 


-05 


54 


58.43 


0.4368 


01 


39 


59.10 


-05 


56 


39.69 


0.4445 


01 


39 


53.30 


-05 


52 


17.88 


0.4486 


01 


39 


59.30 


-05 


53 


30.57 


0.4455 


01 


39 


53.97 


-05 


52 


58.03 


0.4471 


01 


39 


59.51 


-05 


54 


58.52 


0.4522 


01 


39 


56.81 


-05 


54 


24.24 


0.4378 


01 


39 


59.81 


-05 


53 


50.06 


0.4418 


01 


39 


57.27 


-05 


54 


46.19 


0.4311 


01 


39 


59.84 


-05 


54 


31.90 


0.4502 


01 


39 


57.76 


-05 


52 


3.45 


0.4504 


01 


40 


0.13 


-05 


56 


49.62 


0.4379 


01 


39 


59.08 


-05 


54 


15.97 


0.4448 


01 


40 


0.26 


-05 


53 


13.41 


0.4430 


01 


39 


59.98 


-05 


52 


20.03 


0.4468 


01 


40 


0.45 


-05 


53 


52.89 


0.4556 


01 


40 


0.22 


-05 


54 


34.68 


0.4481 


01 


40 


0.48 


-05 


55 


46.69 


0.4464 


01 


40 


0.74 


-05 


55 


1.72 


0.4502 1 


01 


40 


0.83 


-05 


54 


26.79 


0.4412 


01 


40 


0.83 


-05 


55 


2.56 


0.4504 


01 


40 


1.70 


-05 


55 


19.46 


0.4550 


01 


40 


1.10 


-05 


55 


6.69 


0.4405 


01 


40 


2.11 


-05 


54 


50.22 


0.4473 


01 


40 


1.47 


-05 


55 


9.67 


0.8729 3 ' 6 


01 


40 


2.13 


-05 


54 


27.62 


0.4451 


01 


40 


1.96 


-05 


55 


13.94 


0.4542 


01 


40 


2.27 


-05 


56 


38.88 


0.4425 


01 


40 


2.13 


-05 


55 


29.26 


0.8270 b 


01 


40 


2.51 


-05 


55 


10.36 


0.4354 


01 


40 


3.14 


-05 


55 


20.58 


0.45 10 2 


01 


40 


2.75 


-05 


53 


35.74 


0.4405 


01 


40 


3.84 


-05 


55 


21.90 


0.4547 


01 


40 


2.80 


-05 


56 


48.62 


0.4451 


01 


40 


4.39 


-05 


56 


6.87 


0.4398 


01 


40 


2.84 


-05 


53 


49.91 


0.4463 


01 


40 


5.13 


-05 


55 


37.15 


0.4412 


01 


40 


2.86 


-05 


55 


26.48 


0.4398 


01 


40 


6.02 


-05 


55 


4.80 


0.6096 b 


01 


40 


3.24 


-05 


55 


6.47 


0.4508 


01 


40 


6.11 


-05 


56 


27.71 


0.4446 


01 


40 


3.28 


-05 


55 


34.00 


0.4504 


01 


40 


8.40 


-05 


54 


58.29 


0.4231 


01 


40 


3.49 


-05 


55 


22.51 


0.4541 










DEIMOS 




01 


40 


3.89 


-05 


55 


32.22 


0.4514 


01 


39 


53.62 


-05 


55 


6.64 


0.4396 


01 


40 


4.36 


-05 


54 


59.94 


0.4508 


01 


39 


55.39 


-05 


55 


22.26 


0.4541 


01 


40 


4.96 


-05 


56 


10.91 


0.3332 


01 


39 


55.59 


-05 


55 


10.04 


0.4472 


01 


40 


5.11 


-05 


54 


10.70 


0.4434 


01 


39 


55.84 


-05 


54 


12.44 


0.4437 


01 


40 


5.37 


-05 


56 


43.50 


0.4560 


01 


39 


55.93 


-05 


55 


1.28 


0.4497 


01 


40 


6.18 


-05 


54 


11.23 


0.4418 


01 


39 


56.90 


-05 


54 


35.23 


0.9321 b 


01 


40 


8.51 


-05 


54 


55.35 


0.4249 


01 


39 


56.92 


-05 


55 


30.78 


0.4520 


01 


40 


9.21 


-05 


54 


20.65 


0.4417 


01 


39 


57.06 


-05 


55 


22.06 


0.4515 

















Note. — Typical errors in the measured redshifts are 0.0001 to 0.0002. 

1 BCG1, the north-western BCG. 
2 BCG2, the south-eastern BCG. 
3 A 1, the central gravitational arc. 
b Background object. 



onto a mass concentration described by a King profile w ith po c± 
2.6 x 10 25 g cm -3 and r c ~ 210 kpc. |Bradac et al.| ( |2008| ) estimate 
a merger velocity of -2000 km s" 1 for MACS J0025. 4-1 222 by 
scaling the results of Markevitch et al. (2004) to match the velocity 
dispersion and mass of their tamet cluster. In either case, the merger 
velocity is clearly supersonkjj owing to the deep gravitational po- 
tential well. Since the X-ray gas mass of MACS J0140.0-0555 is 
comparable to that of MACS J0025. 4-1 222, and the total masses 
of the two systems are also comparable (discussed below), we can 



4 In the X-ray emitting gas of massive clusters sound travels at a speed of 
approximately 1300 km s _1 . 



expect the merger velocity of MACS J0140.0-0555 to be the order 
of a few 1000 km s -1 . 

The observed small relative radial velocity of the two BCGs 
may be explained by any of the following three scenarios: (1) 
the two clusters are approaching each other along an axis almost 
aligned with our line of sight but are still so far apart that the Hub- 
ble expansion essentially cancels out the aforementioned free-fall 
velocity; (2) the two clusters are colliding along an axis almost 
aligned with the plane of the sky such that the observed difference 
in radial velocity represents only a small fraction of their relative 
velocity in three dimensions; (3) the two clusters are observed near 
turnaround after the primary collision and are close to at rest with 
respect to each other — in this final scenario, the orientation of 
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Figure 7. #ST colour images ((a):1.5' x 1', (b) & (c): 20" x 20"; blue: ACS/F606W, green: WFC3/F110W, red: WFC3/F140W). In (a), the core region of 
MACS J0140.0-0555, grey contours show the adaptively smoothed X-ray emission and red contours show the critical lines at z = 0.87 derived from our mass 
model (see Section [T6) . The BCGs of the two sub-clusters as well as the central gravitational arc are labeled. Redshifts as measured with LRIS and DEIMOS 
are labeled in white in all panels. Strong lensing features compatible with our mass model are marked with blue dashed circles (one in each panel (a), (b) & 
(c)). We only successfully measured the redshift for the feature in (c). The lower left corner in (a) and the upper right corner in (c), as marked by the grey 
dashed lines, were only observed with ACS. The noisy circular region lower left to the center of (c) are dead pixels on WFC3. 




4000 5000 6000 7000 8000 9000 

Wavelength (A) 



Figure 8. Spectrum of Al, the gravitational arc between the two BCGs (see Fig.[7}, as observed with LRIS. The spectra obtained with the blue and red arm 
of the instrument are not flux calibrated and have been median- smoothed to increase S/N. Identified emission lines are labeled. Grey vertical bands mark the 
locations of absorption and emission features in the Earth's atmosphere. We derive a redshift of 0.8729 for Al, which places it well behind MACS J0140.0- 
0555. 



the merger axis is barely constrained at all. Of these three scenar- 
ios, we can firmly rule out the first one, since the required separa- 
tion of several Mpc would render the system a pre-merger. Prior to 
first contact, however, the gaseous cores of the two clusters should 
align well with their BCGs which is in conflict with the observed 
X-ray /optical morphology shown in Fig. 1 . Distinguishing between 
scenarios 2 and 3 is much more difficult based on the data at hand. 
Assuming that turnaround after the primary collision occurs ap- 



proximately at the virial radius (r2oo), the observed separation of 
the BCGs in the plane of the sky of about 250 kpc implies an incli- 
nation of the merger axis of less than 82° under scenario 3 — a very 
weak constraint indeed. Scenario 2 (the one adopted by Bradac et 
|aL] ( [2008l ) for MACS J0025.4-1222) is appealing mainly because 
of its simplicity; however, absent all signs of shocks in the X-ray 
data it is not strongly favoured by the observational evidence. 

Although any conclusions regarding the impact parameter of 
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the collision must necessarily remain somewhat speculative too, 
three pieces of observational evidence support the classification of 
MACS J0140.0-0555 as a post-collision, binary, head-on merger 
by |Mann & Ebeling | < [20T2] >: (1) the excellent alignment of the X- 
ray emission with the line connecting the BCGs of the two subclus- 
ters; (2) the elongation of the X-ray contours along the merger axis 
suggested by that same line; and (3) the absence of two X-ray peaks 
that could be associated with the cores of the individual subclusters 
and the position of the single X-ray peak between the BCGs. 

To summarize, we conclude that MACS J0140.0-0555 is a bi- 
nary merger proceeding at close to zero impact parameter. In view 
of the absence of clear signs of a very recent collision (be it in the 
form of shocks, surviving cluster cores, or pronounced disturbance 
in the ICM), we propose that MACS J0140.0-0555 is observed well 
after the primary collision and quite possibly after turnaround. In 
this scenario the inclination of the merger axis with respect to the 
plane of the sky cannot be constrained beyond stating that it is 
likely to be less than approximately 80°. 

4.2 Cluster scaling relations 

The estimate of the total mass of M(< r 50 o) = (6.8 - 9.1) x 10 14 
M , combined with the global gas temperature of 8.5^q § keV and 
a 0.1-2.4 keV luminosity of L x (< r 50 o) of 6.5 x 10 44 erg s" 1 
(see Section[T2) places MACS J0140.0-0555 within the scattering 
of the Luminosity-Mass and Temperature-Mass scaling relations 
( |Mantz et al.|20T0] ). 

4.3 Implications for the properties of dark matter 

Although the lack of deep HST imaging prevents us from mea- 
suring the mass distribution in MACS J0140.0-0555 using weak- 
lensing techniques, the estimates derived from the measured gas 
mass and from the simple strong-lensing mass model (Sections |3.2| 
and |3.6| > strongly suggest total masses for the two subclusters that 
are compa rable to those of M ACS J0025.4-1222 of « 2 - 3 x 
10 14 M jBradac et al.|2008| >. MACS JO 140. 0-055 5 thus has the 
potential to provide independent constraints on the dark matter self- 
interaction cross section. 



5 SUMMARY 

Our joint X-ray / optical analysis of MACS JO 140.0-0555 based 
on observations with Chandra and HST finds a velocity disper- 
sion of a = 875^q km s _1 (n z = 66) as well as Lx,boi(< 
r 50 o) = 2.1 x 10 45 erg s" 1 , kT - 8.5 keV, and M(<r 50 o) = 
(6.8 - 9.1) x 10 14 M . While the system's X-ray / optical mor- 
phology clearly identifies MACS J0140.0-0555 as a post-collision 
merger, there remains little observational evidence of the colli- 
sion in the form of shock fronts, temperature gradients, or a bi- 
modal redshift distribution of the cluster galaxies. We conclude that 
MACS J0140.0-0555 is a well advanced merger. Although all ob- 
servational evidence is consistent with the merger proceeding at 
minimal impact parameter and close to the plane of the sky, we 
are unable to formally constrain the inclination of the merger axis 
with respect to the plane of the sky, due to the inherent degener- 
acy between (radial) peculiar velocity and Hubble flow. Heuristic 
arguments, however, suggest that MACS J0140.0-0555 is a head- 
on, binary merger along an only mildly inclined axis, observed 
well after the primary collision and possibly after turnaround of the 



merger components. Our strong-lensing analysis of the system con- 
firms the presence of two well separated mass distributions, mak- 
ing MACS J0140.0-0555 a potential candidate for a measurement 
of cr/m, the self-interaction cross section of dark matter. However, 
resolving the two subclusters, which are separated by only 43" (250 
kpc) in projection, would pose a challenge for a weak-lensing anal- 
ysis. 
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